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The pion and proton mass in finite volume 
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We calculate the finite volume effects for the pion and nucleon mass. For the pion mass we present the results 
of a full two-loop calculation in chiral perturbation theory. The outcome shows that the resummed version of 
the Liischer formula we presented in an earlier work does indeed give an excellent approximation to the two-loop 
result. In view of this result we apply the same resummed formula to the nucleon mass. In the nucleon sector 
the extension of the chiral expansion to higher quark masses appears to be more problematic and it is therefore 
more difficult to make reliable predictions for the size of the finite volume effects. We present some preliminary 
numerical estimates. 


1. Introduction 

One of the many useful applications of effective 
field theories is to estimate the size of systematic 
effects in lattice calculations. This concerns both 
the unphysical size of the quark masses as well 
as the finite lattice spacing which, in the frame¬ 
work of the effective field theory, can be seen as a 
controlled modification of the underlying funda¬ 
mental theory. Since one can keep track of this 
modification even at the level of the effective La- 
grangian, one is able to calculate how the physics 
is changed at large distances. A conceptually dif¬ 
ferent problem is to estimate the effects due to 
the finite volume - in this case one is modifying 
the physics at large distances but not touching 
the short-distance physics. Correspondingly, one 
can evaluate the distortion due to the finite vol¬ 
ume by calculating any observable in infinite and 
in finite volume in the framework of the effective 
Lagrangian This evaluation can only be reli¬ 
able if the new scale introduced by the finite vol¬ 
ume can be denoted as a “low energy scale”, i.e. 
if ^ A, where the latter is a typical hadronic 
scale. In order to have a better numerical esti¬ 
mate of how large a volume should be, we further 
observe that the minimum nonzero momentum 
allowed in a box of size L is p = 27r/L, and that 
if one identifies A with (as is usually done 
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in chiral perturation theory) one obtains 

2LF^ » I =k L» Ifm . (1) 

In recent years several analytical calculations of 
finite volume effects within chiral perturbation 
theory (CHPT) have been performed by various 
groups (for a recent review, see [ 2 ]). These re¬ 
sults may be used as a guideline by lattice groups 
when planning their runs and deciding on the size 
of their volumes. Since these effects are gener¬ 
ally small (at the percent level) one could either 
try to estimate the volume for which these ef¬ 
fects are negligible (below one percent, say), or 
alternatively to work in the smallest possible vol¬ 
ume in which the effect can be reliably calculated 
and correct analytically for the finite volume. In 
the latter case it is of course important that the 
uncertainty in the analytical calculation be care¬ 
fully estimated. Moreover, before blindly trust¬ 
ing some theoretical calculations, it is certainly a 
good practice to check them at the claimed level 
of accuracy - before one can save some computer 
time on the volume at least some should be de¬ 
voted to checking the volume dependence. 

In order to estimate the uncertainty in the fi¬ 
nite volume calculations it is necessary to check 
the convergence of the chiral expansion. This re¬ 
quires evaluating at least two terms in the series 
for this specific effect. Since finite volume effects 
start only at the one loop level, a check of the 
convergence of the series requires a two-loop cal- 
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culation. Until recently m no such calculations 
had ever been performed. 

A shortcut which did not require a full two- 
loop calculation but relied on the Liischer formula 
for the masses (and on its extension to decay 
constants [S|), was proposed in m Asymptotic 
formulae d la Liischer express the finite volume 
shift as an integral over an infinite volume ampli¬ 
tude - e.g. the finite volume effect for the pion 
(proton) mass is expressed as an integral over the 
TTTT (ttTV) scattering amplitude. Inserting the tree 
level representation for the latter amplitude in the 
Liischer formula yields the asymptotically domi¬ 
nating term of the one-loop finite volume calcula¬ 
tion. But since the relevant infinite-volume scat¬ 
tering amplitudes are usually known well beyond 
the tree level, one can insert a better represen¬ 
tation than the leading order one and easily go 
beyond the leading order calculation of the finite 
volume effect. 

This has been done in m for the pion mass. 
The somewhat surprising results were: 1. the 
leading order in the chiral expansion receives very 
large corrections at next-to-leading, but reason¬ 
able ones at next-to-next-to leading order; 2. the 
leading exponential term (as given by the Liischer 
formula) is numerically dominating only where 
the finite volume effect becomes negligible. The 
conclusions to be drawn from these results were 
that in order to reliably calculate these effects it 
was necessary to go beyond the leading term both 
in the chiral expansion as well as in the asymp¬ 
totic expansion. The possibility to use the NNLO 
representation for the tttt scattering amplitude 
had however allowed us to show that, apart for 
the anomalously large NLO correction (the rea¬ 
son for this large correction is well understood, 
cf. |H]), the series behaved as expected and for 
not too large pion masses was clearly converg¬ 
ing. The shift from NLO to NNLO was typically 
smaller than the uncertainties in the NLO calcu¬ 
lation due to the low energy constants. 

These results lead us to formulate a simple 
recipe [J] to evaluate finite volume effects beyond 
LO. The recipe can be described as a resummed 
Liischer formula: the integral over the infinite vol¬ 
ume amplitude which appears in the latter is gen¬ 
eralized to one with an index n such that the 


asymptotic behaviour becomes exp(—-ynMTrL), 
and the full result is given by a sum over all n, 
with appropriate coefficients. This recipe yields 
the exact one-loop CHPT result if one inserts the 
LO amplitude in the integral - at NLO the cor¬ 
respondence is not exact anymore, but we argued 
that we expected the resummed asymptotic for¬ 
mula to give the largest part of the full two-loop 
calculation. Moreover we could show that alge¬ 
braically the formula improved the accuracy of 
the asymptotic formula and that corrections were 
of the order 0(e“'^^) with [M = /\/2M,r) 

instead of the with M = of 

the plain Liischer formula. 

Recently, we have completed the full two-loop 
calculation of the pion mass and could confirm ^ 
that the resummed asymptotic formula is an ex¬ 
cellent approximation to the full two-loop results, 
and can therefore be used with confidence also in 
other cases. One interesting application is that of 
the nucleon mass: the latter has been calculated 
to one loop in CHPT m and the result has suc¬ 
cessfully described lattice data, despite the fact 
that these had been obtained for substantial pion 
masses and in rather small volumes. With the re¬ 
summed formula we can go beyond the one-loop 
calculation and better estimate the uncertainties. 
These two new results (which will be published 
soon) will be briefly described here. 


2. The pion mass to two loops 

In Ref. |2I7| we have proposed a resummed ver¬ 
sion of the asymptotic formulae d la Liischer as 
an efficient way to evaluate finite volume effects 
in CHPT beyond leading order. The Liischer for¬ 
mula for the relative finite volume shift for the 
pion mass reads as follows 


= 


M^{L) - 




( 2 ) 




/ OO _ 

dy .F,(i 2 /)e-vWA, 

-OO 


+ 0(e-"^) , 


where J^Tr{v) denotes the infinite volume forward 
(t = 0) TTTT scattering amplitude in Minkowski 
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space, and 


Att — 


( 3 ) 


is the box length in pion mass units. As discussed 
in the introduction, A,r is assumed to be large, and 
the integral in m is only the dominating term in 
an asymptotic expansion in exponentials of mul¬ 
tiples of A^. The estimated error in 0 is de¬ 
termined by M > -\/3 /. In Liischer’s proof 
of the formula the first step is to show that the 
dominating contribution to finite volume effects 
comes from diagrams where only one propagator 
is taken in finite volume. The propagator in finite 
volume is 


GL(a;°,x) = ^Goo(x°,x-hnL) , (4) 

n 


where Goo is the standard, infinite-volume prop¬ 
agator. The term with n = 0 in the sum is sim¬ 
ply the infinite-volume propagator. According to 
the rules for doing CHPT calculations in hnite 
volume one has to proceed as in infinite vol¬ 
ume and simply write all propagators as in 0 . 
As shown by Liischer the dominating term comes 
from graphs where one propagator at a time is 
taken in finite volume. He then showed that in 
this class of diagrams the dominating contribu¬ 
tion comes from the terms with |n| = 1 in the 
sum in 0 . Contributions from all other terms in 
the sum in 0 , however, will also show up in the 
complete CHPT calculation and can be dealt with 
exactly as the first term in the sum. This is the 
origin of the resummation we suggested in m 
and leads to the resummed formula 



where m(n) is the number of integer vectors z 
with z^ = n. It is easy to verify that the term 
with n = 1 is identical to the original Liischer 
formula 0. Notice that although the formula 
contains subleading terms in the asymptotic ex¬ 
pansion, the algebraic accuracy of the formula is 
in principle the same. On the other hand, an 


analysis of the contributions at the two-loop level 
not included in the resummed asymptotic formula 
shows that at this order in the chiral expansion 
0 the formula receives corrections only of order 
exp(—ML) with M = M7r(-\/3-|-1)/-\/2, and so 
is also algebraically improved with respect to the 
original Liischer formula. 

On the basis of this and other arguments we 
have used this resummed formula to make a 
thorough numerical analysis of finite volume ef¬ 
fects for masses and decay constants of the pseu¬ 
doscalar mesons [J]. It is important, however, to 
check the claim that the resummed formula pro¬ 
vides the main contribution of a full NLO calcu¬ 
lation of the finite volume effects, at least in one 
concrete example. We have now done that by 
calculating the pion mass to two loops in CHPT. 

The pion mass in finite volume to two loops can 
be written as 


M^{Lf = M^- Ml = , 

where 

= I^ + I^ + 0{e) , ( 6 ) 


h = 


Ml ^m{n) 


^ IGTr^Air , ^/n 

n—1 ^ 

oo 

J dy L'^(i 2 /) 


L = - 


iMi 


E 


i(n) 


327r3A,r —; \/n 

n—1 ^ 

oo oo /— - — 

n(s+y^)A„ 


Jdy jds 


s + 2iy 


disc[L'^(s,l-biy)] , 


and introduced the abbreviation 

Ml 


( 7 ) 


The expression for E*^^) shows that at the two- 
loop level there are contributions from diagrams 
with only one propagator in finite volume which 
are not included in the resummed formula 0. 
We have denoted these as Ic and shown that these 
can also be represented compactly, as a double in¬ 
tegral over the off-shell tttt scattering amplitude. 
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Figure 1. Relative finite volume effect for the pion 
mass as a function of the pion mass for L = 2, 3 
and 4 fm. 


Contributions from diagrams where more than 
one propagator is taken in finite volume are by 
definition not included in the resummed Liischer 
formula. The whole contribution is a two- 
loop effect which goes beyond the latter and can 
be represented as 

5^(2) ^ 

+ 0{y) , (8) 


3. The nucleon mass with the resummed 
Liischer formula 


In this section we apply the resummed Liischer 
formula to the nucleon mass and first provide the 
relevant analytical expression. The difference to 
the case of the pion mass is that now the ttN 
scattering amplitude, instead of the tttt, has to 
be inserted in the integral. 

The formula for the relative finite volume cor¬ 
rection i?Ar = (TOAr(L)—TOAr)/mAr reads as follows 


34 ^ m{n) 


o 2 -y'n(l-£2)A„ 


47r2 ynK 


(9) 


where 


S’tt — 


2m N 


, D+iy)=mND+{iM^y,0) , (10) 


where the latter is one of the components of the 
elastic ttN scattering amplitude, which is defined 
as follows: 


T(7r“(g)iV(p) ^ 7r-'{q')Nip')) ^ T,,, = , , 

= Sa'aT+ + ^[Ta',Ta]T- . ^ ’ 

Each of the two isospin components, r+ and T~ 
is then broken down into 


=u' 


D^iy, t) - 


AmN 


t) 


( 12 ) 


where gi is the one-loop tadpole graph. The part 
from the sunset type diagrams which can not be 
written in terms of gi is denoted here with A. 
It can be written as double integrals and can be 
evaluated numerically. 

The numerical analysis of this formula is shown 
in Fig. nj It is evident that the contributions 
which go beyond the resummed Liischer formula 
are sizeable only in the region where is not 
large, i.e. where the framework in which we are 
doing our calculations is not reliable anymore. 
The results show that the resummed Liischer for¬ 
mula provides an efficient and reliable way to go 
beyond leading order in the evaluation of finite 
volume effects. 


and each of the amplitudes depends on two kine- 
matical variables: 


t = {q-q')y 

s = {p + qf , 


V = 
U = 


s—u 
4m At 


{p-q'Y 


(13) 


As seen in (Emu, the ttN amplitude is needed 
here in a particular kinematical configuration, 
namely for t = 0 and for v purely imaginary and 
small: the contributions with large values of v 
are suppressed by the exponential weight in the 
integral in ®. It is therefore natural to make a 
Taylor expansion of the amplitude around v = 0 
after having subtracted the pole due to the one- 
nucleon exchange diagram (also called the Born 
term). Such an expansion is in fact already well 
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known in the phenomenology and is referred to as 
the subthreshold expansion. It reads as follows 


D+{u,G) = D+,{v,Q) + D+{u,G) + Dt,{s,u) ,(14) 
where 


i?+(z.,0) 


91n 

2 9 ’ 

rriN 

^00 + *^10^ + ) 


( 15 ) 


and vb = —M^/{2mN). The function D^^{s,u) 
contains the analytically nontrivial part of the 
amplitude. Up to order in the chiral expan¬ 
sion this can be written as a sum of two single 
variable functions: 


Dt,{s,u)=Dtis) + Dt{u) + 0{q^) ( 16 ) 


which admit the following dispersive representa¬ 
tion: 


Dt{s) = 


du' - 


I 


ImD^ (s') 

— u — ie) 


( 17 ) 


where s' = 2m+ m'j^ -I- M^. This representa¬ 
tion shows that, due to the large number of sub¬ 
tractions, the function Df is small near i/ = 0. Its 
contribution to the finite size shift of the nucleon 
mass is negligible. 

This observation leads us to the following ex¬ 
pression for the relative finite volume shift 


Rn = 


3e^ ^ m(n) 

dTT^ ^ y/nX.^ 
n=l 




57rAf^7r^pv('\/^^'7r) ^00^® (V^^tt) 


-I- (i)‘oS^(\/nA^) - dtoB'^{VnX-„) 


where 

-^pv (-^TT ) 


dy- 


-OO 

^OO 


4 + 2/" 


f dy 

J —OO 


( 18 ) 


( 19 ) 


are the relevant finite volume integrals, df^ = 
mNM^'^d^Q and i?Ar,na is the remainder coming 


from the (subtracted) analytically nontrivial part 
of the amplitude. 

If we neglect the contribution RN,na the rep¬ 
resentation dm is very simple and expresses the 
finite volume shift of the nucleon mass in terms of 
only a handful of physical observables: the pion 
and proton masses, tuat, the pion-nucleon 
coupling constant and the three subthreshold 
parameters df^ (the latter are not directly observ¬ 
able, but can be obtained from data with some 
theoretical treatment, cf. CH) If one knows the 
low-energy constants (LEG) which appear in the 
chiral representation of these quantities, one can 
predict their quark mass dependence and there¬ 
fore the finite volume shift i ?at as a function of 
the quark masses. 

An explicit representation for the quark mass 
dependence of the quantities which appear in 
Eq. (tUHl up to 0{p*) can be found in Un¬ 
fortunately, however, our knowledge of the LEG 
which appear in there is much worse than that 
for the LEG of tttt scattering. In short we can say 
that the quark mass dependence pf the nN scat¬ 
tering amplitude, even in the subthreshold region, 
i.e. far away from the physical singularities of 
the scattering amplitude, is not very well known. 
We stress that for the physical value of the quark 
masses we do have reliable phenomenological in¬ 
formation about the quantities which appear in 
Eq. Ijl 811 . and that the finite volume effects can be 
therefore evaluated with rather small uncertain¬ 
ties, cf. m The problem is the extrapolation to 
higher quark masses. 

We have investigated this question and adopted 
the following strategy: we fit the physical values 
of the five quantities g a and the three sub¬ 
threshold parameters dio (cf. ^J) as well as the 
lattice calculations of (cf. cuni) and gA na 
for Mtt ~ 0.5 GeV. This gives us seven constraints 
and allows us to fix seven of the LEG appearing in 
the formulae for the mass dependence. On some 
of the other LEG there is phenomenological infor¬ 
mation coming from ttN or ttN tt ttN scatter¬ 
ing - whenever possible we use this also. Some 
of the LEG, however, remain unconstrained and 
have to be estimated on a purely theoretical basis. 
The details of this analysis, including a discussion 
of the numerical values of the various LEG, will 
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be discussed elsewhere M- 

An important point we wish to emphasize here 
is that our analysis indicates that the Born term is 
dominating even at higher values of pion masses. 
This means that the most important information 
we need, concerns the quark mass dependence of 
mN and For the former we rely on available 
lattice calculations (although at somewhat high 
pion masses), whereas for the latter we rely on 
the Goldberger-Treiman relation: 

= ( 20 ) 


and fix the LEG dis from the phenomenology. 
This gives a mild quark mass dependence in the 
relation - much more important is the quark mass 
dependence of gA and uin- Also for gA we rely 
on available lattice calculations (again at some¬ 
what high pion masses), which indicate a mild de¬ 
pendence of the axial charge on the quark mass, 
and try to interpolate the physical and the lat¬ 
tice value with the chiral representation. In this 
case, however, the phenomenological information 
about the LEG (which determine the quark mass 
dependence of gA) would actually allow us to pre¬ 
dict the latter m 


9 a — 9 + ( 4di6 ~ 


5" 


167r2E2 J 

(c3 - 2cA)gM^ 




{l + g^)gM^ 


SirmF'^ 


67rE2 


0{M^) 


( 21 ) 


The constant dig has been determined from 
ttN —> ttttN scattering measurements [HI to be 
between —1 and —2 GeV“^, with about 1 GeV“^ 
uncertainty^, whereas the combination C 3 — 2 c 4 
can be determined from the subthreshold parame¬ 
ters of ttN scattering m and comes out to be also 
negative, with a value around —9 GeV“2. Unfor¬ 
tunately the prediction indicates a very strong 
quark mass dependence already at very small 
quark masses and makes it difficult to extrapolate 
to the region where the lattice data are obtained: 
the O(M^) correction is large and negative and 


2 The published values of die have later been corrected in 
the final version of the PhD thesis of Nadia Fettes (see also 
the discussion in El)' The numerical change, however, 
does not change the qualitative picture. 


the 0{M^) large and positive. Of course one can 
choose the value of gA in the chiral limit such 
that the value for the physical pion mass comes 
out right, but as soon as one goes to higher pion 
masses, gA tends to explode. The lattice data, on 
the other hand, indicate that at ~ 0.5 GeV, 
gA is somewhat lower, but not by much, than the 
physical value. We have to conclude that either 
the phenomenological determination of the LEG 
appearing in gA is unreliable, or that the region 
of quark masses where the chiral expansion works 
for gA is very small. 

The determination of dig is indeed subject to 
large uncertainties (of the order of 50 to 80% ac¬ 
cording to the estimates of Fettes, which do not 
seem to include any systematic effects), but those 
of C 3 and C 4 are more solid and indicate a small 
region of convergence of the chiral expansion for 
gA ■ The presence of these large corrections has in¬ 
deed been known since a long time (cf. EH) 
a possible solution of the discrepancy with the 
observed mass dependence (or the lack thereof) 
in lattice calculations has already been proposed 
by Hemmert, Procura and Weise CHI In this 
paper the A resonance is included explicitly in 
the calculation and it is shown that if one stops 
at order in the small-scale expansion (SSE) one 
obtains a relatively steep dependence only close to 
the chiral limit, whereas at higher quark masses 
gA is a flat function of M^. While this analysis 
identihes a mechanism that provides a change of 
behaviour between the very low and the middle 
pion mass region, it does not provide a fully sat¬ 
isfactory understanding of the mass dependence 
of gA, in our opinion. The constant dig occurs 
also in the SSE as an independent LEG, in fact 
unrelated to the A resonance. In m the value 
of dig is also taken from the work of Fettes et al. 
m this means that even in the SSE the term 
of order is a large correction already at the 
physical pion mass, and that the size of this cor¬ 
rection grows fast with the pion mass. The con¬ 
tributions from the A, which are also large, com¬ 
pensate these large corrections and give a rather 
flat behaviour up to 0.5 GeV. However, since the 
A and the dig contributions are physically un¬ 
related (also in the framework of the SSE), this 
compensation appears to be the result of a fine 
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Figure 2. Relative finite volume effect for the 
nucleon mass as a function of the pion mass. 


tuning between the value of die and the A reso¬ 
nance couplings. The analysis in m shows that 
the A contributions could tame the strong depen¬ 
dence of gA on the quark mass, but does not yet 
explain why this happens. 

This situation makes the numerical study of fi¬ 
nite volume effects for the nucleon mass on the ba¬ 
sis of Eq. m problematic. We proceed nonethe¬ 
less and rely on a simple interpolation formula 
between the physical and the lattice value (at 
around M^r ~ 0.5 GeV) of gA - for all other quan¬ 
tities we rely on the chiral expansion. Our results 
are shown in Fig. [21 and indicate that although 
with large uncertainties, the finite volume cor¬ 
rections for the nucleon mass can be calculated 
with this method. The figure contains two curves 
(with the corresponding uncertainties): the dot- 
dashed curve indicates what one obtains if one 
stops the chiral expansion of the ttN scattering 
amplitude at 0{p^), whereas the solid line gives 
the result obtained from the full 0{p^) scattering 
amplitude. In the latter case in the amplitude 
there are more LEG to be determined and cor¬ 
relations among them are generated by the fits, 
but the net result is that the uncertainties become 
larger than if one stops at 0(p^). The comparison 
clearly shows that the direct statistical estimate 
of the uncertainty at 0{p^) is quite optimistic. 


We conclude by comparing to the analysis in 
cni. The algebraic relation between the two ap¬ 
proaches is fully understood: if we expand the 
resummed Liischer formula Eq. llTKll to NLO we 
reproduce exactly the result of Ref. cni. The par¬ 
tial inclusion of higher orders obtained if we do 
not expand Eq. CHI but insert the formulae of 
the 0{p‘^) ttN scattering amplitude, on the other 
hand, allows us to check the convergence of the 
series. Fig. (21 shows a reasonable behaviour for 
L = 2 fm and Mt^ < 0.5 GeV. In Ref. the 
comparison between lattice data and the GHPT 
calculation in finite volume has been performed 
even for smaller volumes and higher pion masses 
and has shown a surprisingly good agreement. 
The significance of this agreement can however 
only be assessed if one has a good estimate of the 
uncertainties, something which had not been at¬ 
tempted in m- We can do this now and compare 
our results to the same lattice data, which were 
obtained for ttin at ~ 0.55 GeV and volumes 
L <2 fm - a region where we would not trust our 
formulae from the start. The comparison is shown 
in Fig. 121 where one can see that if one stops at 
NLO one goes indeed through the data, in agree¬ 
ment with cni. The inclusion of higher orders, 
however, spoils the good agreement, because for 
the central values of the LEG we have obtained 
from our fit we come closer to the LO curve. The 
band which gives our estimate of the uncertainty 
however shows that neither the agreement nor the 
disagreement are of any significance, because the 
uncertainties at these high quark masses and for 
such small volumes are simply too large. 

4. Conclusions 

We have discusssed the results of a full two- 
loop calculation of the pion mass in finite volume 
in GHPT We have shown that contributions 
which are not included in the resummed asymp¬ 
totic formula we proposed some time ago m are 
tiny and can be neglected in the region of pion 
masses and volumes where GHPT can be trusted. 
This check provides further evidence to our claim 
that the resummed asymptotic formulae are the 
most efficient way to evaluate these finite volume 
effects and in view of this there appears to be lit- 








Figure 3. Finite volume effect for the nucleon 
mass as a function of L for Mjr = 0.545 GeV. The 
lattice data are from Ref. d The shaded band 
represents the uncertainty of the calculation. 


tie need to try and improve the analysis in [Jj for 
masses and decay constants of the pseudoscalar 
mesons. 

The resummed formula can also be applied 
to non pseudo-Goldstone bosons, and a particu¬ 
larly interesting application concerns the nucleon 
mass. This has been calculated in finite volume 
to one loop in GHPT cni and compared to lattice 
data, although these had been obtained for rather 
heavy pion masses and in small volumes. Some¬ 
what surprisingly, the GHPT calculation success¬ 
fully described the volume dependence of these 
data m- We have described here a numerical 
analysis based on the resummed asymptotic for¬ 
mula and shown how higher order corrections be¬ 
come more and more sizeable as the pion mass 
increases. With the higher order corrections, also 
the uncertainties in the calculation increase. In¬ 
deed we have shown that in the region of pion 
masses and volumes where data are available the 
uncertainties are larger than the effect itself and 
concluded that the observed agreement between 
lattice data and the GHPT NLO calculation is 
accidental and of little significance. 
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